High-Field Shubnikov-de Haas Oscillations in the Topological Insulator Bi2Te2Se. 
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We report measurements of the surface Shubnikov de Haas oscillations (SdH) on crystals of the 
topological insulator Bi2Te2Se. In crystals with large bulk resistivity (~4 f^cm at 4 K), we observe 
~15 surface SdH oscillations (to the n — 1 Landau Level) in magnetic fields B up to 45 Tesla. 
Extrapolating to the limit 1/B — s- 0, we confirm the ^-shift expected from a Dirac spectrum. The 
results are consistent with a very small surface Lande (;-factor. 
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I. INTRODUCTION 

In Topological Insulators, the surface electrons occupy 
helical Dirac states in which the spin is locked perpen- 
dicular to the momentum [H-Q]- In three-dimensional 
examples, the topological surface state was observed h 
angle-resolved photoemission spectroscopy (ARPES) 
[1]. Scanning tunneling microscopy (STM) has also been 
applied extensively j9l4ll|. In transport experiments, 
quantum oscillations of the surface electrons have been 
observed in BizTeg [H, and in (Bi,Sb)Se3 [H. The 
Quantum Hall Effect was also observed in a thick film of 
strained HgTe ,14]. However, in the Bi-based materials, 
progress has been slowed by the small surface conduc- 
tance C relative to the bulk term G^. We report mea- 
surements on crystals of Bi2Te2Se in which C/G^ ^^1 
and SdH oscillations with large amplitudes are observed 
at high fields. By tracking the Landau Level (LL) ex- 
trema towards the quantum limit, we observe directly 
the i-shift that distinguishes the Dirac spectrum from 
the Schrodinger case. Our results address the question 
whether the spin-Zeeman energy affects the LL sequence 
in the quantum limit. 

Landau quantization of the surface Dirac spectrum was 
previously observed in Bi2Se3 by STM Nonethe- 
less, high-i? transport experiments to approach the quan- 
tum limit are important to search for novel states. In ad- 
dition, accurate determination of the i-shift associated 
with the Berry phase provides the best test for whether 
the SdH oscillations arise from surface topological states 
or bulk states (this requires a large B to reach the n = 
1 LL). 

In a magnetic field B normal to the surface, the Dirac 
states are quantized into Landau Levels (LLs). As B is 
increased, sequential emptying of the LLs leads to oscilla- 
tions in . We follow the customary practice of defining 
the "index field" Bn as the field at which the Fermi en- 
ergy Ep lies between two LLs, i.e. at the minima in G^ 
(see Sec. [TTI- A plot of the integers n vs. l/Bn gives a 
nominally straight line with slope equal to the FS cross- 
section Sp. 

Our interest is in the limit l/Bn ^ 0. In the 
Schrodinger case, there are n filled LLs below Ep when 
the field equals i?„ (as defined). By contrast, in the Dirac 



case, we have n+ ^ filled LLs between Ep and the Dirac 
point (at E = 0). The important additional ^ arises be- 
cause the conduction band and the valence band each 
contributes half of the states that make up the n — LL. 
Hence, as 1/_B ^ 0, the plot of vs. n intercepts the 
n-axis at the value 7 = for the Dirac case, whereas the 
intercept 7 = (mod 1) in the Schrodinger case. The ^- 
shift was experimentally verified for the Dirac spectrum 
in graphene, and expressed equivalently as a Berry-phase 
TT-shift il5|]. 



II. RESISTIVITY MAXIMA OR MINIMA? 

The index field Bn clearly plays the key role in pin- 
ning down the shift in the index plot. Here we wish 
to discuss the question of determining Bn when surface 
and bulk carriers co-exist [16]. In the bismuth-based 
systems (and other 3D topological insulators), the two- 
dimensional electron gas (2DEG) on the surface is in in- 
timate contact with bulk electrons which conduct a sig- 
nificant fraction of the applied current. By contrast, the 
entire current is carried by the 2DEG in graphene and 
GaAs heterostructures. When Ep falls between adjacent 
LLs in the QHE regime of graphene, both the 2D con- 
ductance Gs and resistance Rxx attain a deep minimum 
(this follows from Ry^ ^ Rxx)- 

However, when a large, parallel bulk conduction chan- 
nel exists (the case here), the observed conductance ma- 
trix is the sum 



G, 



7 1 



(1) 



where G?, is the bulk conductance matrix. As the mo- 
bility of the bulk carriers /ifc is very low (50 cm^/Vs), 
bulk SdH oscillations are not observable even at 45 
T. The additivity of the conductances in Eq. [T] im- 
plies that the index fields still correspond to minima 
in Gxx- However, because the bulk G^^ is dominant, 
the observed resistance now attains maxima at Bn (i.e. 
Rxx = Gxx/[Glx + Gly] - 1/G,,). We find that it is 
least confusing to work with Gij because its components 
are additive. The results reported here provide an exper- 
imental verification of this point. 
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FIG. 1: (color online) A representative set of the observed 
resistivity p and Hall coefficient Rh vs. T in Bi2Te2Se (sam- 
ples identified by the numbers). The magnitudes of p and 
Rh at 4 K vary considerably between annealed samples. In 
all samples, the carriers are predominantly n-type at 4 K. The 
change in sign of Rh near 56 K reflects the thermal activation 
of bulk hole carriers across a gap of 50 mV. The largest SdH 
amplitudes are observed in samples with p > 4 Qcm at 4 K. 
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In many experiments, however, the Hall response is not 
available. One may still use the SdH oscillations in the 
resistance Rxx, provided i?„ is identified with its maxima. 
If the wrong choice is made (identifying i?„ with minima 
in Rxx), a spurious intercept will appear for carriers 
with a Schrodinger dispersion. 

A second issue we address is the strength of the Zee- 
man energy. Strict particle-hole symmetry implies that 
it is unshifted in energy. On the other hand, a large 
Zeeman energy q^bB may lead to high-field distortion 
of the SdH period {g is the surface Lande g-factor and 
/is the Bohr magneton). The in- field STM experiments 
[Tol [Tl| have shown that the n — Q LL is unshifted up 
to 11 Tesla. This test can be extended to much larger 
B in transport experiments, but early SdH experiments 
had limited resolution [l^, [11] ■ Values of g as large as 
76 have been inferred from low-field SdH oscillations in 
BiaTeaSe [13 • 



III. EXPERIMENTAL DETAILS 

The large density of Se vacancies (electron donors) in 
BiaSea leads to an n-type semi-metal with a sizeable car- 
rier density (rif, ~ 10^^ cm^'^). By contrast, as-grown 
crystals of BiaTes are p-type because of Te-Bi exchange 
defects. In the hybrid material BiaTcaSe, the Se ions oc- 
cupy the innermost layer in each quintuplet layer. This 
appears to suppress both vacancy formation and Te-Bi 
exchange defects. Two groups have found that surface 
SdH oscillations are observed in n-type crystals with 
greatly reduced nf, [l^ [3 ■ Details of the crystal growth 



FIG. 2: (color online) The resistance (per square) Rxx and 
Hall conductance Gxy in Bi2Te2Se (Sample 4). Panel (a) 
shows the SdH oscillations in Rxx vs. B at fields above 11 T 
at T = 0.7 and 2.5 K. At 40 T, the peak-to-peak amplitude is 
17 % of the observed resistance. The Hall conductance Gxy 
at 0.7 K is plotted in Panel (b). In both panels, the envelope 
is the smooth curve passing through the extrema points. The 
background curve (dashed curve) is determined as the average 
of the envelope curves. 

for our samples appear in Ref. f23|. 

Even in carefully annealed crystals, large variations 
in the values of n?, and the observed resistivity p are 
found [2(|. Figure [T] shows traces of p vs. T for a repre- 
sentative set (Samples 1, 2 and 3). At 4 K, p varies from 
1 to 6 ricm. Although all these samples exhibit SdH os- 
cillations, the amplitudes are largest when p > 4 ilcm at 
4 K. 

As shown, the Hall coefficient Rh changes from p to 
n-type as T decreases near 56 K. We have found [21]] that 
the Hall behavior results from the thermal activation of 
holes into the bulk valence band across a "transport" gap 
At ^ 50 mV. Previously, we showed [l^ that the surface 
conductance ff* in BiaTeaSe involves carriers with a high 
mobility ps of 2,800 cm^/Vs, whereas the residual bulk 
conductance (from an impurity band) involves n-type 
carriers with much smaller mobility (/if, ~ 50 cm^/Vs). 
The magnitudes of inferred from kp and ps confirm 
that the SdH oscillations are from surface states. Ando's 
group has shown in field-tilt experiments that the SdH 
period is consistent with surface states [3]. Helical sur- 
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face states in an isolated Dirac band have been observed 
by spin-resolved ARPES 

The large variation in p may be understood by esti- 
mating the number defects. If we assume that each de- 
fect (either Se vacancies or Te-Bi exchanges) contributes 
a carrier, the observed rif, (3x10^^ cm~'^ in Samples 1 
and 2) corresponds to a defect density of a few parts 
in 10^ [i^l- This stringent constraint implies that fluc- 
tuations at this level lead to pronounced variations in 
rih and p. Even in optimally annealed crystals, separate 
portions of an exposed surface can display different p- 
T profiles. In addition, aging of the surface results in a 
gradual decrease in the amplitude of the surface quantum 
oscillations with time (roughly by a factor of 2 over a few 
weeks for crystals sealed in Ar atmostphere and stored 
in dry ice) . These factors are problematical for high- field 
transport experiments. 

To improve the odds, we cleaved crystals ~30 minutes 
before loading the high-field cryostat. Each crystal was 
contacted by 3 pairs of leads so that both the resistance 
tensor i?^ can be measured over distinct segments. Be- 
cause the 45-Tesla field cannot be reversed, we employed 
the reciprocity technique of Ref. 23] to extract both R^x 
and Ryx- 
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IV. QUANTUM OSCILLATIONS 

We report measurements to fields of 45 T in Samples 
1 and 4 (in which Rh = -137 and -52 cm^/C, respec- 
tively, at 4 K). The large, well-resolved SdH oscillations 
in these samples provide an opportunity to investigate 
the specific issues in the quantum limit. As shown in 
Fig. the peak-to-peak SdH amplitude in the resis- 
tance Rxx in Sample 4 grows with B until it accounts 
for ~17% of the total resistance. Because conductances 
are additive, it is expedient to convert Rij to the conduc- 
tance Gxx = Rxx /[Rxx + Ryx] ^-nd the Hall conductance 
Gxy = Ryx/[Rlx + RlxV Gxy is plotted in Fig. Hd. 
Using the envelope of the oscillations (faint curves) , we 
locate the midpoint between adjacent extrema to define 
the background. 

After removing the background, we isolate the oscilla- 
tory components l^Gxx and l^Gxy which we plot versus 
l/B in Fig. [3] The conductance l^Gxx and Hall con- 
ductance AGxy are plotted in Panels (a) and (b), respec- 
tively (both normalized to the quantum of conductance 
e^/h). The fit of the oscillations (see Sec. E]) yields a 
surface mobility of 3,200 cm^/Vs and a metallicity pa- 
rameter kpi = 30. The interesting phase shift apparent 
at low B is discussed later. 

FigureU^ plots the minima of AGxx versus n (solid cir- 
cles) . In addition, the maxima of A.Gxx have been plotted 
as open circles (shifted by i). The best-fit straight line 
gives a Fermi cross-section area Sp of 48.5 T. A similar 
plot based on the extrema of the Hall conductance AGxy 
is shown in Fig. [Jb- The minima in —AGxy correspond 
to n-|-|;, since the derivative —dAGxy/dB has minima at 



FIG. 3: (color online) The oscillatory component of the con- 
ductance AGxx (Panel a) and the Hall conductance AGxy 
(Panel b) in Sample 4 plotted against 1/B (T = 0.7 K). The 
two quantities are normalized to /h. The fit of the oscil- 
lations (see Fig. ^ yields a surface mobility of 3,200±300 
cm^/Vs, with kpi ^ 30. In Sample 4, G" accounts for ~ 19% 
of the total conductance at 4 K. Note the phase shift at low 
B. The LL indices n = 1,2,3 are indicated for the minima of 
AGxx- 

n. The value of Sp found from AGxy (47.3 T) is consis- 
tent with the previous value within our resolution. The 
values of n = 1,2,3 at the minima of AGxx are noted in 

Fig. m^- 

In order to fix the intercept 7, we expand the scale in 
Fig. 21;. The best-fit straight line passing through the 
six extrema of AGxx intercepts the n-axis at the value 7 
= -0.61±0.03. Similarly, the high-field extrema of AGxy 
are plotted in Fig. The intercept for the best-fit 

line occurs at 7 = -0.37±0.03. Within our uncertainties, 
these intercepts are significantly closer to the ideal value 
7 = — i than or 1. Hence, the high- field results provide 
transport evidence for a Dirac spectrum for the surface 
states. 

Although we do not observe quantized Hall steps in 
Fig. (the oscillatory component rides on a large tilted 
background contribution from the bulk Hall current), it 
is interesting that the peak-to-peak amplitude swing of 
AGxy is ~0.8 e'^/h per surface for n — 1, which is of the 
order of the quantized Hall conductance value. 

In Sample 1, the amplitudes of the observed SdH os- 
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FIG. 4: (color online) The index plots of 1/-B„ vs. the integers 
n in Sample 4. In Panel (a), Bn is obtained from the minima 
of AG XX- In Panel (b), the index field B'^ is inferred from the 
minima of —AGxy B'„ is plotted against n+ j, where the j 
shift arises because the minima in dAGxy/dB align with the 
minima in AGxx- We expand the scale in Panels (c) and (d) 
to show the intercepts more clearly. In Panel (c), the solid 
straight line is the best fit to the extrema fields for n <3. The 
dashed line is the best fit to all the extrema field shown in 
Panel (a). The sketch shows Ef in relation to the filled LLs 
(solid color) in the Dirac spectrum when B — 42.0 T (arrow). 
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FIG. 5: (color online) The oscillatory component AGxx vs. 
1/B (Panel a) and the index plot of 1/Bn vs. n (Panel b) in 
Sample 1. The intercept 7 of the best-fit line is -0.45±0.02. 



dilations are considerably weaker (Fig. [5^). The index 
plot of l/Bn vs. n fits a straight line that intercepts the 
n-axis at 7 = -0.45±0.02, again consistent with a Dirac 
spectrum. 

The expanded plot shows why intense fields are needed 
to fix 7 reliably. By accessing the n — ^ index at 45 T 
(Figs. |lt,d), we have reduced considerably the "spread" 
of intercepts caused by the measurement uncertainties: 
an intercept 7 = may be safely excluded. A more 
subtle point is the slight curvature of the index plot. In 
Fig. Ht, if we extrapolate the best-fit line (dashed) using 
the total data set from 3 to 45 T, its intercept yields 
-0.78, nearly exactly between -1 and -5. By contrast, 
the best-fit line (bold) to the high-field extrema for n <3 
yields an intercept (-0.61) closer to -i. This implies that 
the index curve 1/Bn vs. n develops a slight curvature 
in intense fields. (The curvature accounts for the low-_B 
phase shift apparent in the single-frequency fit in Fig. 

I3i.) 

A possible cause of curvature is the spin-Zeeman en- 
ergy. When that is included, the Hamiltonian is 

H = VFh-(TXTv-^B-(T (2) 

where n the unit vector normal to the surface, cr are the 
spin Pauli matrices, and tt = p — eA is the momentum 



p of the electron in a vector potential A. The LL energy 
is given by 



En = ±^2nhvleB + {g^iBB/2)^. (3) 

The energy of the n = LL increases linearly with B 
instead of being unshifted. For a large g, the plot of 
1/Bn vs. n will deviate from a straight line as 1/B — > 

0. In our experiment, we have tracked the LLs to rt = 

1. The weak deviation from a straight line in Fig. 

is inconsistent with values of g substantially larger than 

2. More importantly, however, the observed deviation is 
opposite in sign to that predicted by Eq. [31 As we do not 
see evidence for a deviation caused by a large g-factor, 
we conclude that the the g factor of the surface states 
in Bi2Te2Se are not significantly greater than 2 in the 
quantum limit. 

V. SURFACE CARRIER MOBILITY 

In general, it is very difficult to separate G" from 
reliably even at B = 0. Shubnikov de-Haas (SdH) os- 
cillations - when measured with sufficient resolution - 
provide a powerful way to tease out the surface conduc- 
tance. Analysis of the SdH amplitude vs. B yields the 
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estimate that, in Sample 4 at B=0, G" accounts for 
19% of the total observed conductance. These values are 
similar to those obtained in an earlier sample, which had 
a slightly larger kp (0.047 A'^) 

The mobility provides a strong, quantitative argument 
that the SdH oscillations arise from surface states. Sup- 
pose for the sake of argument that the oscillations arise 
from bulk states. The SdH period is then to be identi- 
fied with a 3D Fermi sphere of radius kp — 0.038 A~^, 
or a 3D carrier density of 1.86x10^^ cm^"^. With this 
density, the inferred mobility gives a 3D resistivity p\, ^ 
1.1 mJlcm at 4 K. Instead we measure p to be 5 J7cm. 
The large discrepancy (factor of 4,500) firmly precludes 
a bulk origin for the SdH oscillations. 



FIG. 6: (Color online) The oscillatory component of the con- 
ductance ^Gxx in Sample 4 at 0.7 K (solid curve) and the fit 
to Eq. [5] using only one frequency (dashed curve). 



scattering rate and the surface mobility /ij (equivalently 
the mean- free-path €). Also, the period of the oscillations 
yields hp. With /i^ and kp known, we then obtain the 
zcro-i? value of = G"* using 



= {e^/h)kp(.. 



(4) 



To focus on the SdH oscillations, we first determine the 
envelope curves passing through the extrema of the oscil- 
lations as explained in Fig. 2 of the main text. The oscil- 
latory component /S.Gxx is obtained by subtracting from 
Gxx the background, defined as the curve lying between 
the envelope curves. (We remark that AG^x does not 
account for all of the surface conductance. By construc- 
tion, its field-averaged value {AGxx)b vanishes. Hence 
we must have AGxx < Gxx-) 

To fit the oscillatory component AGxx, we employed 
the standard Lifshitz-Kosevich expression (2^ 
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sinh A 



2itEt. 
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, (5) 



with A = 2t:'^ k bT /fvjjc and A^ — 2Tr'^kBTp)/hujc, where 
ujc is the cyclotron frequency and the Dingle tempera- 
ture is given by To = h/ [iirkBT), with r the lifetime. 
For 2D systems, we may write the SdH frequency as 
2ttEpB /{huJc), which simplifies to A-K^hns/e, with the 2D 
carrier density rig = fc|,/47r (per spin). Equation [5] may 
be employed in a Dirac system if we write the cyclotron 
mass as rac = E/vp. 

As shown in Fig. [6l we obtain a reasonably close fit 
to the observed oscillations (bold curve) using just one 
frequency. The optimal fit yields for the 3 adjustable 
parameters the values kp = 0.038 A~^, ip = 0.657r and 
Tp) = 8.5±1.5 K, which implies a surface mean-free-path 
i = 79±8 nm and mobility ^i, = e£/hkp = 3,200±300 
cm^/Vs. The metallicity parameter kp£ equals 30. We 



VI. CONCLUSIONS 

The Dirac-like topological surface states detected in 
ARPES and STM experiments present a host of new op- 
portunities for transport experiments especially in high 
magnetic fields. In bulk crystals, the presence of bulk 
carriers complicate transport studies. As shown here, 
quantum oscillations provide a powerful way to isolate 
the surface carriers and to determine their mobility and 
kp£. The index plot of the integers n versus l/Bn can 
be used to confirm the 7r-shift associated with the Berry 
phase of the surface electrons, which leads to an inter- 
cept -i in the limit 1/B — > 0. To access LLs at n = 1 (or 
lower), we have employed fields up to 45 T. The results in 
Figs. 21 and [5] provide direct confirmation of the existence 
of the intercept expected from a Dirac disperion. 

The resolution attained here provides experimental 
verification of the point that the intercept is observed 
only when i?„ is identified with minima in Gxx or max- 
ima in Rxx- (For contrast, we note a recent report (25| in 
which a -i intercept was obtained in high-B measure- 
ments on exfoliated crystals of Bi2Te3. However, be- 
cause Bn was inferred from minima in the resistivity, 
it seems that the -i intercept actually implies a Berry 
phase that is zero, consistent with SdH oscillations from 
conventional bulk carriers. ) 

The linearity of the index plot in Figs. |4] and [5] show 
that the Lande (/-factor is small {g ~2). The n = LL is 
unshifted even at 45 T, consistent with STM experiments 
taken at 11 T [l3,[iii. 

Finally, we comment on the results in the large- B limit. 
In Fig. the last maximum in AGxx (at i? ~ 40 T) 
corresponds to n = ^ (see arrow in the index plot in Fig. 
I3J:). At this field, the Fermi energy Ep is aligned with 
the center of the n = 1 LL, as sketched in the inset in 
Fig. lit. In our indexing scheme, there is 1 filled LL 
between Ep and the Dirac Point, with ^ of the filled 
states from the unshifted LL at the Dirac Point). Hence, 
these results provide rather firm evidence for this ^ -shift 
in the limit 1/B — > 0. As the inset in Fig. HJ; implies, the 
interesting states in the n — LL in Sample 4 become 
experimentally accessible in fields higher than 45 T. 
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